1 



On 
OS 
t— i 



> 
o 

oo 

o 

On 



X 



Monte Carlo computation of the effective potential for the three 
dimensional Ising system 

Jae -Kwon Kim and D. P. Landau 

Center for Simulational Physics, The University of Georgia, Athens, GA 30602 

Using a novel finite size scaling Monte Carlo technique, we calculate the four, six and eight point renormalized 
coupling constants defined at zero momentum for the three dimensional Ising system. Our values of the six and 
eight point coupling constants are significantly different from those obtained from other methods. 



1. INTRODUCTION 

Effective action is a useful formalism to de- 
scribe long distance, low momentum behavior 
of a system. In this work we are interested in 
the Monte Carlo calculation of the the effective 
potential for the three dimensional (3D) Ising 
model. The Ising model is equivalent to the sin- 
gle component X4> 4 theory in the limit A — > oo, so 
that the use of the standard weak coupling per- 
turbation method is plagued because of the strong 
coupling nature of the model. Hence a variety of 
new analytical methods have recently been de- 
veloped (Tj-Q. Due to the non-perturbative na- 
ture of the problem, Monte Cairo methods have 
also been used. However, standard Monte Carlo 
measurements of relevant physical quantities suf- 
fer from high statistical noise(|]. In this work 
we employ a novel finite size scaling (FSS) tech- 
nique combined with a cluster flipping Monte 
Carlo algorithm in order to calculate the coef- 
ficients of the Taylor expansion of the effective 
potential, which are closely related to the dimen- 
sionless N-point renormalized coupling constants 
(RCC). 

In terms of the renormalized field ip R the effec- 
tive potential can be written as 
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where m R and represent respectively the 

renormalized mass (inverse correlation length (£)) 
and the N- point RCC defined at zero momen- 



tum. It is easy to show thatQ 
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Here Z v represents the field strength renormal- 
ized factor (Z v = X m R j with \ denoting mag- 
netic susceptibility), and Wjv is the Fourier trans- 
formed Appoint connected Green function of zero 
momentum. 

2. THE MONTE CARLO METHOD 

Since we are concerned with the low momen- 
tum properties of the model when it is in the 
scaling regime (but not at criticality), our goal 
is to calculate the thermodynamic values (the in- 
finite volume limit value) of the RCC measured 
at zero momentum in the limit m R — > 0. ffjj is 
expressed in terms of the combinations of the var- 
ious connected Green functions, so that the sta- 
tistical noise in its numerical estimation becomes 
increasingly large as N increases. For a fixed N, 
we empirically observe that the statistical error 
increases drastically for approximately L/£ > 4. 

To illustrate we measured various physical 
quantities at an arbitrary (inverse) temperature 
with increasing linear size of the lattice L. Ta- 
ble(l) clearly shows that the statistical error in 
the measured values of g R increases with N for 
any given value of L. For each quantity, the cor- 
responding value increases monotonically with L 
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Table 1 Table 2 

Size dependence of the various physical quantities Size dependence of the various physical quantities 
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1285(10) 
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18.3(26.4) 


72 


19.3(1) 


20.3(4) 


1515(85) 


16.3(3.3) 
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5.62(1) 


26.0(4) 


2089(278) 















(see also Table(2)) until it becomes L indepen- 
dent, and can be regarded as the correspond- 
ing thermodynamic value. The standard theory 
of FSS shows that the thermodynamic limit is 
reached beyond a certain universal value of L/£, 
regardless of the value of the temperature. 

The correlation length in Table(l) measured 
under the condition > 4 seems to be almost 
L independent. However, this is certainly not 
the case for other renormalized quantities which 
continue to increase significantly until — 6. 

Note also that for the six and eight point renor- 
malized couplings the relative statistical errors in- 
crease drastically for larger values of L, whereas 
those for the correlation length remains negligi- 
bly small. Evidently, the requirement of larger 
L for the traditional Monte Carlo measurement 
of the thermodynamic values of the RCC makes 
it virtually impossible to measure these quanti- 
ties precisely in a deep scaling regime due to the 
large statistical fluctuations. 

We would like to stress that the number of our 
Monte Carlo measurements is very large. For 
each L we typically generated tens of millions of 
the single cluster configurations. For L=28, for 
example, Monte Carlo measurements were taken 
every 20th configuration out of a total of 2 x 10 s 
configurations that were generated; yet, the rela- 
tive statistical errors for the g^p and gjjr amount 
to approximately 4% and 25% respectively. 

In order to overcome the difficulty closer 
to T c , we make use of a new FSS function 



(Qa(x(L, i)))@], defined by the expression 

A L (t) = A(t)Q A (x(L, t)), x(L, t) = a(t)/L. (1) 

Here A^it) represents the quantity A measured 
on a finite lattice of linear size L at a reduced 
temperature t, with its corresponding thermody- 
namic value A(t) scaling as A(t) ~ t~ p . 

The technique is especially useful for our pur- 
pose, because it enables us to extract accurate 
thermodynamic values based on the Monte Carlo 
measurements on much smaller lattices. We just 
outline the FSS extrapolation technique used in 
this work. For a detailed explanation, we refer 
the readers to Ref. 

1) . For a certain to, measure j4l(£o) and 
x(L,to) = ^L(to)/L for increasing L. 

2) . Determine A(to) by measuring j4l(£o) such 
that it is L independent. 

3) . Fit (x(L,t ),A L (t )/A(t )) data to a suitable 
functional form. In this work we used the ansatz, 
Q(x) = 1 + Cl e- l /' x + c 2 e- 2 / x + .... 

4) . For any other t, choose a suitable L, measure 
the value of x(L, t) = (,l/L and A^it), and inter- 
polate Q(x(L, t)). 

5) . Extract A{t) by plugging A L (t) and Q(x(L, t)) 
into Eq.(l). 

Our t is (3 = 0.217 and 0.220, and our largest 
value of L for the extraction of the thermody- 
namic values up to (3 = 0.2213 is 72. 

3. RESULT AND DISCUSSION 

At (3 — 0.217 we have been unable to determine 
the precise value of the thermodynamic value of 
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Table 3 

Thermodynamic values of the four and six point RCC extracted by the FSS technique for some temper- 
atures over 0.217 < (3 < 0.2213. 

— 8 0217 0\219 (1220 0.2206 0.2210 0.2212 0.2213 

g^ 25.6(4) 25.3(1.1) 24.7(5) 24.6(7) 24.5(3) 24.5(3) 24.4(3) 

g ( R ] 2100(100) 2074(144) 1982(157) 2006(157) 1949(146) 1966(135) 1952(127) 



g R due to the large statistical noise in the mea- 
surements for L=32 and 36, but we estimate that 
g<g\/3 = 0.217) ~ 2100(100). With this value, 
we have been able to determine the thermody- 
namic values of the g R up to (3 = 0.2213 (where 
£ ~ 29.2(1)). Of course, all other quantities in- 
cluding correlation length have been determined 
without any ambiguity ||] since their thermody- 
namic values can be accurately measured. We 
summarize our results for the RCC in Table(3), 
from /3 = 0.217 to f3 = 0.2213. 

We observe that g^ R decreases quite slowly as 
£ increases, but eventually the value becomes sta- 
blized around 24.5(2). g^ R has a tendency to de- 
crease as well, but within the statistical fluctua- 
tions it remains a constant for (3 > 0.2210, i.e., 
g^ R ~ 1956(145). The statistics of our g R data 
at j3 = 0.217 are not good enough for the ap- 
plication of our FSS technique. At (3 = 0.2210, 
however, we have been able to measure g R quite 
precisely (Table(2)) by varying up to 3.73, 

where we find g R = 1.63(33) x 10 5 . Again we ob- 
serve a systematic increase of g R ^ with increasing 

L, and estimate g^ ~ 1.75(25) x 10 5 . Our results 
indicate that all the RCC "scale" with the cor- 
responding values of the critical exponents zero. 
Since the susceptibility and the correlation length 
already "scale" over this regime of the tempera- 
ture ||, such pathological behavior as a sudden 
change of them closer to criticality is very un- 
likely. 

(4) 

Our value of g R is in reasonable agreement 
with other estimates, but our estimate of g^ R is 
more than 40% larger than those obtained from 
the measurement of the probability density of 
the order parameter 0, and from average action 
method and linked cluster expansion Q. The 
disagreement in g^ R is more conspicuous; our 



value is approximately three times larger than 
those estimated from the other methods. The 
Monte Carlo measurements for the probability 
distribution in Ref. Q, however, were made un- 
der the condition ~ 4, which is probably 
too small. Our result contradicts the claims in 
the recent literature that the effective potential 
can be very well approximated by the terms up to 
ip 6 . On the other hand, our estimates of the 
renormalizcd couplings are significantly smaller 
than those estimated by the dimensional expan- 
sion||. It has been noticed lately, however, that 
the higher order terms in the expansion become 
negative p. 
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